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ABSTRACT . 


The electromagnetic dispersion relation for two counterstreaming 
ion beams of arbitrary relative strength flowing parallel to a dc magnetic 
field is derived. The beams flow through a stationary electron background 
and the dispersion relation in the fluid approximation is unaffected by the 
electron thermal pressure. Magnetic effects on the ion beams are 
included but the electrons are treated as a magnetized fluid, m^-^0. 

The dispersion relation is solved with a zero net current condition applied 
and the regions of instability in the k-U space (U is the relative velocity 
between the two ion beams) are presented. These results are extensions 
of Kovner's analysis for weak beams. The parameters are then chosen 
to be applicable for parallel shocks. We find that unstable waves with 
zero group velocity in the shock frame can exist near the leading edge 

V. 

of the shock for upstream Alfven Mach number's" greater than -5—5. - -It- i-s 

suggested that this mechanism could generate sufficient turbulence within 
the shock layer to scatter the incoming ions and create the required 
dissipation for intermediate strength shocks. 
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INTRODUCTION 


The instabilities that result from the presence of two inter- 
penetrating ion streams immersed in an electron background may 
provide the mechanism for collisionless momentum coupling between 
the two streams. The presence of a weak magnetic field can have a 
significant effect on the physical process. First, gradients in the 

magnetic field give rise to electron currents that can drive ion acoustic 

1 2 

waves unstable and increase the effective collision frequency; ’ 

3- 6 

second, when the propagation is perpendicular to the magnetic field, 
the magnetic field can inhibit the electrons from shorting out ion plasma 
oscillations for wavelengths long compared to the electron gyroradius; 
and third, whistler mode-ion beam interactions are likely to be 
important, and the existence of whistler modes depends on the presence 
of a magnetic field. It is this third aspect, and in particular, its 
application to collisionless oblique and parallel shock wave structures, 

that is of interest in" this- paper. - - _ 

Significant among recent advances in magnetic streaming in- 
stabilities and their applications to collisionless oblique shock wave 

7 8 

structures are the so-called *'beam- cyclotron'' (BC) ' and "beam- 
9 

whistler” (BW) instabilities. These are characterized by possible 
interaction between one of the two covmter streaming ion acoustic beam 
modes and the electron cyclotron or whistler mode (Fig. 1), While 
the BC instability does result in electron heating, it nevertheless maybe 





inadequate as a mechanism for ion- ion momentum coupling, for this 

instability can be stabilized by electron heating or by resonance broadening 

7 

or by ion trapping. On the other hand, BW instabilities do result in 
turbulence which can suitably couple the two ion beams. The application 
of BW interactions to collisionless shock structures is especially 
attractive since they are found to be operative over a substantially 
larger range of Mach numbers than are the so-called "magnetic counter- 
streaming ion-ion"^ ^ and "modified"^' ^ two stream instabilities 
associated with perpendicular shocks. 

It is interesting to note that for parallel propagation (k || B^) 
the ion acoustic modes of the counter streaming ions uncouple from the 
whistler mode in the linear BW dispersion relation cited by Biskamp 

Q 

and Welter. Thus, in their assumed frequency range 
(flci and being the ion and electron cyclotron frequencies), the 

whistler modes are stable. On the other hand, if one relaxes the 
frequency restriction to lower frequencies, i,e. , 0<|cj| « then, 

even in the case of parallel propagation, a beam mode can, in fact, 
enter into unstable interaction with the whistler leading to coupling of 
the two ion beams. The object of this paper is to analyze this instability 
especially for parameters that can be related to parallel shocks. 

In Section II we present an analysis of the whistler mode 
instability arising in a collisionless plasma configuration consisting of 
two ion beams counter streaming along the magnetic field in the rest frame 
of the background electrons. In Section III we shall determine when 
these unstable whistler modes can grow to large amplitude inside a 
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shock front propagating along the magnetic field. These calculations 
are carried out in the shock frame at different upstream Mach numbers 
and at different stations within the shock thickness. The result of these 
calculations is a map showing the range of Mach numbers as a function 
of position into the shock for which unstable whistlers can stand inside 
the shock. A summary of the results is presented in Section IV. 



LINEAR ANALYSIS 


We consider an unbounded collisionless plasma consisting of two 
ion beams of densities n^^ and n 2 counterstreaming with velocities 
and parallel to a steady magnetic field B taken to lie along the positive 
z-axis. We assume the background electrons to be at rest. The charge 
neutrality and zero current conditions give 

n, + n- = n (1) 

1 Z e 

and 

n, V- + n^ V. - 0 (2) 

1 ^1 2 ~2 ' 

where n^ is the number density of electrons. For the low frequency case 
that we are interested in, it is crucial that the zero current condition (2) 
be met for arbitrary values of the beam strength^^ (0 :< n^/n < 1). ' The 
electron motions are modeled by fluid equations which are valid for wave- 
lengths long compared with the electron gyroradius while the ions are 

assumed to be cold.* The dielectric tensor-for-suc-h a plasma -is -easily 

calculated and the appropriate electromagnetic streaming dispersion rela- 
tion for plane waves exp [i (k. r - cut)] for the case of parallel propagation 
(k II B^) takes the form 

where 
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In the above, CO = (47rn e /m.) ' and CO = (47rn e /m ) ' are the s th 

s 1 ^ ^ 

ion beam (s = 1,2) and electron plasma frequencies respectively. We note 

that within the fluid approximation, the dispersion relation (3) is unaffected 

by electron thermal effects for the case of parallel propagation. 

We now make the low frequency approximations that the displacement 

current is negligibly small (CO « kc) and that the electrons have zero mass 

( |coj « ^ g)- Introducing the dimensionless beam density r\ = n^/n^ 

2 2 

and CO . = Z CO , the dispersion relation (3) simplifies to 

pi ps 


22_ (l-ri)(C0 - k. V ) Ti(a)-k.V) 

^ 1_ ^ 


= 0 


(4a, b) 


CO 


pi 




Cl 


CO - k- V, + a . 

^ ^1 — Cl 


CO - k* + Q . 
^ ~2 — Cl 


where Eq. (4a) refers to the upper signs and Eq. (4b) refers to the lower 
signs. In what follows, Eq. (4b) corresponding to the lower sign will be 
dropped from further consideration since its solution can be recovered from 
that of Eq. (4a) by the substitutions CO ^ - CO and k -k. It is convenient 
^at.this-point to introduce the relative drift velocity U '= Vp“"so'that wd 
have V, = - nU and = (l-n)U where use has been made of the zero current 
condition (2). With these substitutions Eq. (4a) simplifies to 


- 6 - 



U) 


= 0 (5) 
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Now, for 0 < T| <0.5, it is clear that the ion stream 1 has more 
the role of ’'plasma’* while Z has more the role of "beam. " When 
0. 5 < T| <1, these roles are interchanged. However, as the interchanges 

T| (1 - T|) , k- U - k* U 


leave Eq. (5) unchanged, the analysis of the range 0. 5 < rj <1.0 would be 
redundant. It therefore suffices to investigate the range 0 < q < 0. 5. 

For Tj = 0, we recover the usual whistler (Wjf) and ion-cyclotron dispersion 
branches given analytically by 


/2^ci i 

z 1/2 

where C. = (B /47rm.n ) ' is the Alfven speed. For weak beams 
A o 1 e 

(ti « 1), the intersection points P, Q, R (see Fig. 2) between the ion beam 
mode 


CO = ( 1 -Ti) k* U - fi . 


(7) 


and the and IC branches of (6) suggest possible strong interactions near 
these points. Indeed Kovner's analysis" shows that the beam-IC interaction 
is stable while the beam - interactions are unstable with growth rates 


Im 


.1/2 


(k^ 


- 

Cl 


U = p. Q), 
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D5618 


Superposition of the Whistler (W_j_) and Ion Cyclotron (IC^) 
Modes for -q = 0 With a Weak Beam Mode. The intersection 
points R, P, and Q represent regions where strong interactions 
can occur. 
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where, from (6) and (7), kp and depend on the slope (= U/C^) of the 
beam line. We note that the beam-cyclotron (BC) intersection shown in 
Fig. 1 does not appear in Fig. 2 due to our approximation m^ ^ 0. 

We have extended Kovner's weak beam study by solving (5) for 
arbitrary values of and drift Mach number By expanding 

the dispersion relation (5) for small k, one readily finds 


Im Cl) 


k [n (1 - -n) m/ 



SO that these long wavelength modes are unstable only if 


> 

d 


't 1 (1 - Tl)'^ 


( 8 ) 


This condition is thus sufficient for instability, but is not necessary as we 

shall see below. We should emphasize that due to the application of the 

condition (2), we are dealing only with zero-current instabilities and thus 

the stronger, current-driven instabilities often treated are excluded. 

The dispersion curves and growth spectrum from (5) are plotted 

in Figs. 3a, b, c and d for values of r| = 0. 1 and 0. 5 and - 3 and 7. 

In plotting these * -kU- has -been -taken-to be negative with _a view to the 

application in the shock problem discussed in the following section. We 

find that for r| = 0. 1, the peak growth rate Cl)" (Cl) = CO' 1 iCO") increases 

only slightly from Cl)"~ 0. 32 fi . at = 3 to Ct)" 0. 38^2 . at Mj =7. For 

T| = 0. 5, the peak growth rate increases more markedly with rising from 

CO" ~ 0. 34S2 . at M, = 3 to Cl)" 0. 72fl . at M. = 7. For n = 0. 5, the disper- 

ci d ci d ' ^ 

sion curves are symmetric in k as it is impossible to distinguish the beam 
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Fig. 3 The Dispersion Relation (5) for: a) r\ = 0. 1, - 3, 

b) p = 0. 1, Md = 7; c) p = 0. 5,_Md ^0. 5, M d = 7. 

^Note“that'the'Wrfi^ar'scale~should be multiplied by l/’5 

in order to obtain the correct value for Im OJ. 
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f 


from the plasma while for t| = 0, 1, the beam-like mode can be clearly 
identified. We also note that in Fig. 3a the region near k = 0 is stable in 
agreement with (8) but that there is instability for a finite range of negative k. 

The regions of instability in k - U space for r| = 0. 1 and 0. 5, 
obtained by setting the discriminant of the cubic equation (5) equal to zero, are 
shown in Figs, 4a and 4b respectively. Also shown by dashed lines is the 
generalization for arbitrary t| of the Kovner lines. These lines give the 
values of k and U (for a given r\) that simultaneously satisfy (6) and (7) 
and correspond to the intersection points P and Q in Fig. 2. It is easy 
to show that the Kovner lines always lie in the region of instability and that 
they have a minimum at pj = (4/3) 1. 15 for which 2. 6/( l-r|). 

This value corresponds to the relative drift such that the beam line in Fig. 2 
is tangent to the branch, i. e. , the points P and Q coalesce. 

For large drift (M^ °o) there are unstable modes for k in the range 


2 kC 

“d - -- „ - i 72 : h I I H <*-’'> “d 

' ) Cl 




The width of this unstable spectrum is thus 





(i--n)u 


Q . 

Cl 


(9) 


and the peak growth rate is approximately given by CO" — ri ' ^ci' 

These expressions are in agreement with the numerical results shown in 
Figs. 3 and 4. 


- 11 - 





APPLICATION TO SHOCK WAVE STRUCTURE 

The origin of the magnetic field turbulence associated with colli- 
sionless shock waves has not been definitely established. If the instability 
discussed in the preceding section has relevance for parallel shock wave 
structure, then the unstable waves must be able to grow to large amplitude 
inside the shock front. This requirement suggests that the phase or group 
velocity of the whistler in the shock frame should be either zero or very 
small in addition to the wave being unstable. Wave packets with zero group 
velocity can be expected to remain in the shock front for a sufficiently long 
time to reach large amplitude and generate turbulence. In the present 
work, we consider only parallel shock waves for which the sound speed 
Cso = rV/c(T^ -f Tj^)/m^)^^^ and the Alfven speed C^q in the unshocked gas 
are equal. The external magnetic field in this case remains constant and 
parallel to the direction of plasma flow across the shock layer. We thus 
exclude any consideration of switch-on" shocks that occur for 

^SO ^ ^AO- - _ . . _ 

We adopt the Mott-Smith formalism which supposes that the shock 

layer is already a two- streaming environment for the ions. Therefore the 
ion distribution function in the rest frame of the shock front is assumed 
to be a linear combination of the upstream (u) and downstream (d) velocity 
distribution functions, i. e. , 

f. (z, v) = n^(z) 6(v - VJ + n^(z) 6(v - V^). (10) 
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For simplicity, we take to decrease linearly from the upstream 

value at the leading edge of the shock to zero at the trailing edge while 
we assume (z) to increase linearly from zero at the leading edge to 
its downstream value at the trailing edge. A schematic of the density 
and velocity distributions is sketched in Fig. 5. The electrons, on the 
other hand, are treated as a single fluid moving with an average drift * 
velocity, (z), determined by the zero current condition 


V 


e 


(z) 


(z) n. (z) 

■ , , V + ^ , V- 

n (z) n (z) 

e e 


( 11 ) 


= [1 - 11 (z)] + Tl(z) 

where the fraction of the total density represented by the downstream 
species r| (z) is now a function of position through the shock. 

The appropriate dispersion relation is easily obtained from (5) by 
doppler shifting the frequency co by the amount k. (z) which is a 
function of the position in the shock. The instability properties are 
unchanged and depend only on the relative drift velocity U = 


(which is independent of the position in the shock) and r\ (z) (which is 

dependent on the position in the shock). It is convenient, at this stage, 

to specify all the properties and criteria in terms of the upstream 

Alfven Mach niimber the fraction ^ of the distance 

through the shock. From (1), we have 

n (z) = (1 - C) n + C n , (12) 

e u°o ^ doo ! L 

where ^d°o I’epresent respectively the uptstream and downstream 

densities at infinity. Thus 
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( 13 ) 


L = 

^ T1 + T (1 - Tl) 

where T is the compression ratio across the shock which can be obtained 
from the Rankine-Hugoniot relations (with CgQ = C^^q) as 


4M 


T = 


(Y + 1) 


(Y - 1) M^+ 2 M^+ 3 


(for Y = 5/3). 


(14) 


The local Alfven velocity decreases through the shock as the plasma 
density increases and is given by 

1/2 


= '=A0 


3 


M^ + 3 + 3?, (M^ - 1) 


while the local drift Mach number M t is 

d 




V - V 
u d 


3(M^ - 1) 


4M 


A 


1 + 


3^ (M^ - 1) 


M^+3 


(15) 


( 16 ) 


With the help of (11) to (16), the dispersion relation in the shock 
frame can be written as 

0)^ + a (k) + b (k) Gl) + c (k) = 0 

s s s 


where 


CO = CO + k . V (z) 
s ^ ~ e 


and a, b, and c, are known functions of k, and The group velocity 

of the wave 9 co’ / 9k where co' is the real part of co is given by 


8co 

9c0g / 9k = - 


3k + 


2 

i 

3k 


9a ^ , 9 (a + b) , 9 (a b -c) 


3k 


24o)' ^ + 16 a Cxj' + 2(a^ + b) 
s s 


( 17 ) 


for the unstable modes. 
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The unstable regions in k - M^space are as before easily obtained 
by setting the discriminant of (17) equal to zero. These are presented in 
Figs. 6a and 6b for two different values of These figures are different 

from Figs. 4a and 4b since the value of r\ varies for different values of 
for constant We see that near the front of the shock there is a 

wide range of unstable wave numbers for M^> 3. The growth rate is 
of the order ^ However, most of these waves have positive group 
velocities and are swept downstream with the flow before they have 
had a chance to grow significantly. Only unstable waves with nearly 
zero group velocity can be expected to grow to large amplitude. In 
Fig. 6a, the values of k for zero group velocity waves in the unstable 
region are shown by a dashed line. Thus, Mach numbers in the range 
6. 8<M^< 9. 1 should be able to produce turbulence at ^ > 0. 05 according 
to the linearized analysis presented here. Such waves do not exist for the 
parameter ranges shown in Fig. 6b at ^ >0. 2. 

Figure 7 shows the range of upstream Alfven Mach numbers for 

which unstable whistlers will stand in the shock as a function of distance 

into the shock. Two immediate conclusions can be made. The first is’ 

that the mechanism can be important only for M . > 5. 5. The second is 

that these waves are confined to the leading edge of the shock; unstable 

waves near the middle or trailing edge will be swept downstream. This 

13 

result is in agreement with laboratory experiments that measured 
magnetic field fluctuations in lower Mach number shock waves. 
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STABLE 



at a) ?, = 0. 05 and b) ^ 0. 2. The dashed line corresponds 

to values of k for which the group velocity of unstable 
waves is zero in the shock front. 






SUMMARY 


We have derived the electromagnetic dispersion relation for 
two counter streaming ion beams in a stationary electron background. 

The dispersion relation which is unaffected by electron pressure terms 
in the fluid approximation has been solved numerically to find the regions 
of instability in k - M^ space. It is found that for beam- whistler 
modes to be unstable, it is necessary that the relative drift Mach 
number be greater than 2,0 although this condition is not sufficient for 
instability. The bandwidth of k over which the unstable modes exist 
is shown to decrease inversely with the relative drift velocity for 
large relative drift velocities. 

By modeling the shock structure in a Mott-Smith formalism 
as consisting of streaming upstream and downstream ions, the linear 
analysis was applied to the study of shock wave structure propagating 
along the magnetic field. Invoking the criterion that only unstable waves 
with zero group velocity will have time to grow within the shock 
thickness, we found the range of Mach numbers for which some unstable 
beam- whistler modes could grow to large amplitudes. We suggest that 
this may be a possible mechanism to generate turbulence that many 
experiments have shown to be associated with shock wave structures. 

The plasma turbulence in turn can provide the explanation for anomalous 
dissipation. However, in our model, the unstable beam- whistler modes 
can stand near the leading edge of the shock only for upstream Mach 
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numbers greater than 5.5. 

The most serious limitation of the present analysis is the restriction 

k II B . For waves propagating at an angle to the magnetic field, the 
<~«w< o 

dispersion relation becomes much more complicated and electron thermal 
effects have a significant interaction with the magnetic modes treated here. 
The analysis presented here ought to be extended to include the stability 
properties of these obliquely propagating waves. The assumption of cold 
ions ought to be relaxed particularly since one of the ion species in this 
study has been taken to be the shocked ions. Ion thermal effects could 
have a significant effect on the stability properties of counter- streaming 
ion beams. The results of the stability analysis for hot streaming 
plasmas will be reported later. As a final comment, it should be 
emphasized that this instability depends on having included the effect 
of the magnetic field on the ions. If in (4), the instability disappears. 
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